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Abstract 

This paper studies the optimal bit allocation for shape-gain vector quantization of wireless channels 
in multiuser (MU) multiple-input multiple-output (MIMO) downlink systems based on linear precoding. 
Our design minimizes the mean squared-error between the original and quantized channels through 
optimal bit allocation across shape (direction) and gain (magnitude) for a fixed feedback overhead per 
user. This is shown to significantly reduce the quantization error, which in turn, decreases the MU 
interference. This paper makes three main contributions: first, we focus on channel gain quantization 
and derive the quantization distortion, based on a Euclidean distance measure, corresponding to singular 
values of a MIMO channel. Second, we show that the Euclidean distance-based distortion of a unit norm 

_ 2B S 

complex channel, due to shape quantization, is proportional to 2 J «-i, where, B s is the number of 
shape quantization bits and M is the number of transmit antennas. Finally, we show that for channels 
in complex space and allowing for a large feedback overhead, the number of direction quantization bits 
should be approximately (2M — 1) times the number of channel magnitude quantization bits. 

Index Terms 

MIMO Broadcast Channels, Limited Feedback of CSI, Optimal Bit Allocation, Sum Mean Squared 
Error, Shape-Gain Product Quantization. 

I. INTRODUCTION 

The availability of channel state information (CSI) at the transmitter significantly improves 
the performance of multiuser (MU) multiple-input multiple-output (MIMO) systems [[0 — [J3J . 
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Specifically, CSI is essential for effective communications in the MU downlink. In frequency 
division duplex systems, in order to provide the base station (BS) with CSI, the receivers need 
to quantize the CSI and feed the quantized information back to the BS. Clearly this feedback is 
an overhead to the system and, therefore, must be limited to an acceptable level. 

This paper focuses on limited-feedback MU MEMO systems, where a single BS communicates 
with multiple receivers and each user can potentially receive multiple data streams. Specifically, 
we restrict our analysis to systems using linear precoding [4J. Our main goal is to present an 
efficient quantization scheme for these systems. For this purpose, we use sum mean squared error 
across all data streams as the design objective and focus on quantization issues by assuming 
perfect channel estimation at the receiver (user) side and a noiseless, delay-free, feedback link 
to the BS. 

This work is mainly motivated by the fact that the performance of limited-feedback MU-MIMO 
systems is very sensitive to the quality of the CSI available at the BS. Without accurate CSI, the 
quantization error common performance measures saturate in the high signal-to-noise ratio (SNR) 
regime because the BS cannot completely pre-cancel the multi-user interference 0, flU, & It 
is therefore essential to design a limited feedback system such that the CSI quantization error 
is minimized. The optimal design of channel quantization in MU-MIMO systems is, therefore, 
the core objective of this paper. 

Most of the works in limited feedback literature focus on either channel magnitude quantiza- 
tion lEO-lHl or direction quantization 0], ED, OH, JED but not both. However, in MDVIO systems, 
optimizing the precoder at the transmitter depends on both the channel magnitude (also known 
as the channel gain) and the phase of individual channel entries (the channel direction or shape). 
The authors of lfT2l specifically showed that one needs both channel gain and quantized direction 
information to achieve multi-user diversity gain. However, the gain information was assumed 
perfect in [[T2|. In general, joint vector quantization (VQ) of channel magnitude and direction 
is a very complex task [13]. To reduce the design complexity, the authors of lfl4l investigate 
independent quantization of the channel gain and shape and develop optimal bit allocation across 
gain and shape of real channel vectors using spherical codes. The authors of lfT5l . lfT6ll also use 
such a product codebook and solve for the optimal bit allocation to minimize the average transmit 
power with quality-of- service constraints. Such a structure has several practical advantages and 
provides an analytically tractable framework to optimize the limited feedback [[T4|. [fT5l . ifTTl . 
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We adopt a similar approach where channel gain and shape are independently quantized. 

The optimal quantizer depends on the transmission scheme and performance measure used. 
Due to its simplicity and efficiency, we adopt an eigen-based combining (EBC) approach to 
precode the data [Q]|, CO- We study quantization of CSI to minimize the sum mean squared 
error (SMSE) over all data streams received, a popular measure in the MU MIMO downlink flU, 
[fT8l . As shown in our earlier work [19], there is a one-to-one relationship between the SMSE 
objective and the variance of the quantization error. The current work, therefore, focuses on 
optimizing the bit allocation across shape and gain given a budget for feedback overhead per 
user. Once this bit allocation is optimized, one can use our earlier work in [fT9l for designing the 
limited-feedback system. To the best of our knowledge, the problem of optimal bit allocation 
in shape-gain vector quantization to minimize the SMSE of a multiuser MIMO system has not 
been investigated before. 

This paper makes two key contributions: 

1) We show that the quantization distortion of a uniformly distributed unit-norm vector in 
C M is upper bounded by: K s x 2~ 2M - 1 , where M is the total number of transmit antennas, 
B s is the number of shape quantization bits and K s is a constant that does not depend on 

B s . 

2) We also show that, for channels in complex space, the optimal number of channel direction 
quantization bits should be approximately (2M — 1) times the optimal number of channel 
magnitude quantization bits. 

Numerical simulations suggest that the proposed bit allocation laws provide a substantial 
improvement over full shape quantization or full gain quantization in the SMSE and bit error 
rate (BER) performance of a multiuser MIMO system. 

Notation: Lower case letters denote scalar values while lower case bold face letters represent 
column vectors. Upper case boldface letters denote matrices. The superscripts and {-) H 
denote the transpose and conjugate transpose operators respectively. tr[-] denotes the trace 
operator. I is reserved for the identity matrix whereas 1 represents the column vector with all 
one entries. diag(ai, ■ ■ ■ , a n ) denotes the diagonal matrix with diagonal elements ai, • • • , a n ; 
whereas diag(A 1 , ■ ■ ■ , A n ) represents the block diagonal matrix with the matrices A 1; • ■ ■ , A n 
on its main diagonal. || • ||i denotes the L\ norm of the vector. E(-) and S(-) denote statistical 
expectation and surface area respectively. 
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The remainder of the paper is organized as follows. Section [TT] outlines the limited feedback 
MIMO system model and the corresponding shape-gain product VQ structure. Section Hill derives 
the distortion measures and provides the optimal bit allocation solution; in general, proofs are 
deferred to appendices. This section also presents the linear precoding algorithm that incorporates 
the optimal bit allocation policy. Section [IV] presents results of numerical simulations illustrating 
the theory developed. The paper wraps up with some conclusions in Section |V] 

II. System Model 

We begin by developing the system model for linearly-precoded MU-MIMO system followed 
by the model for CSI feedback and the product shape-gain quantization structure. 

A. MU MIMO System Model 

Consider a single base station equipped with M transmit antennas communicating with K 
independent users. User k has N k antennas and receives L k data streams. All data streams are 
independent of each other. Let L = J2k Lk, N = J2k Nk- To ensure resolvability, we require 

L < M and L k < N k . 

Let U G C MxL denote the global precoder, the columns of which are unit-norm. Similarly, 
let P G lZ LxL denote the diagonal power matrix whose entries are the powers allocated to 
individual streams. Let, P max be the total available power; we require tr[P] < P max . The data 
vector x = [xx, ....,xl] T = [xf,x|\ . . . ,x^-] T , includes all L data streams to the K users. The 
N k x M block fading channel, H^, between the BS and user k is assumed to be flat. The global 
channel matrix is H^, with H = [Hi, ...,H fc ]. The elements of channel entries are assumed to 
be zero mean complex Gaussian random variables with unit variance. User k receives 

y£ L = HfUv / Px + n fc , (1) 

where n k represents the zero mean additive white Gaussian noise at the receiver. User k, in 
order to estimate its own transmitted symbols from y k L , forms 

x fc = A fc Vf yf L , (2) 

where V^A^ is the N k x L k decoder matrix for user k. The columns of Y k G C NkXLk are unit 
norm while A k = diag ( X kl , A& 2 , • • • , \ kL J G TZ LkXLk contains the gain variables that normalize 



the received data. Although the gain variables at the receiver side do not affect the signal-to- 
interference-plus-noise ratio (SINR), they play an important role in the error performance of 
transmissions that include amplitude modulation, e.g., quadrature amplitude modulated systems. 
Figure Q] illustrates the proposed downlink system. 

Let VA be the N x L block diagonal global decoder matrix, V = diag (Vi, \r) 6 C NxL 
and A = diag (A 1 , A K ) 6 TZ LxL . Overall, 

x = AV^Uv^x + V^n 

= F^WPx + VV (3) 

where, n = [nf , nj, . . . , n^-] T . For the ease of representation, we define the M x L matrix 
F = HVA with F = [f 1; . . . , i L ], The vectors fi, . . . , are the effective Mxl vector downlink 
channels of the individual data streams. 

The MSE of the i th data stream of the k th user is given b)Q, 



e k,i ~ & 



) (x k)i - x k}i ) H ^ . (4) 



The min-SMSE optimization problem is: 

K L k 

min V" V"efef; subject to tr[P] < P max , ||u £ || = ||v £ || = 1, (5) 

p.ILV.A z — ' z — ' 

' k=l i=l 

To solve this problem, it is computationally efficient to use a virtual dual uplink [4J. In this 
uplink the transmit powers are Q = diag [qi, ■■, qi] for the L data streams, while the matrices 
U and V remain the same as before. 



B. Feedback Model: 

As mentioned earlier, we use an eigen-mode strategy Q, 0. According to this strategy, the 
k th user estimates its own channel Hfc and uses a set of dominant singular values and singular 
vectors of H fc as A fc and V fc respectively. 

Since the user is aware of H^, and A&, it can form the product matrix F& = H^V^A^, 
whose columns act as the effective vector downlink channels for the data streams. Each user 
quantizes its effective vector downlink channel based on an Euclidean distance measure and feeds 

'Note that we, interchangeably, index streams as being the £ th of L streams overall or the i th stream of the k th user. Any 
one-to-one mapping between the two notations is acceptable. 
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back the quantized channel to the BS. Details of the CSI quantization policy will be described in 
the next section. To model the effect of quantization, we consider the following relation between 
the original and the quantized variables, 

fk,i=%,i + fk,i or F = F + F. (6) 

Here, denotes the effective vector downlink channel of the i th stream of the k th user. F 
comprises L effective channel vectors with the original channel directions and channel gains. F 
denotes the L quantized feedback vectors. The matrix F represents the quantization error. 

The BS assumes that the quantization error matrix F has M x L independent identically 
Gaussian distributed (i.i.d.) elements with zero mean and a variance of a 2 E /M, where a% is the 
quantization error variance associated with each quantized vector and is defined as, 



al = E 



(7) 



If — f II 

I l k,i l k,i 1 1 

By using the optimal P and U, the minimum SMSE takes the following form [[191 : 

(8) 



where, 



SMSE = L-M+(a 2 + ^pn**) tr [J" 1 ] , 

J = FQF H + (a* + ^pmax) Im- (9) 

where Q is the virtual uplink power allocation matrix. 

Equations in ® and © show that the SMSE is directly related to the quantization error o\. 
The limited feedback system design problem can therefore be formulated as minimization of the 
quantization error variance subject to a fixed feedback overhead. 



C. Shape-Gain Product Quantization Model 

We intend to find the optimal bit allocation for quantizing the effective vector downlink channel 
ffc j. From now on, we will use z to represent the effective vector downlink channel to simplify 
the notation. According to the eigen-based receiver structure assumed in this work, z represents 
the product of a singular value of the channel matrix and its corresponding singular vector. 

Let z be the quantized effective vector downlink channel and let C = {c 1; c 2 , • • • , c Ntot } denote 
the codebook of quantized channels. Here, N tot = 2 B are the total number of quantization levels 
using a total of B bits. This codebook is simplified to a product codebook. Fig. [2] illustrates the 



7 



product codebook operation based on independent quantization of gain and shape. Let, z = gs 
where, g = ||z|| 2 and the unit-norm s = z/||z|| 2 denote the gain and shape of the channel 
respectively. The BS is provided with the quantized information z = gs, where g and s denote 
the quantized gain and shape respectively. 

Let B g and B s denote the number of bits allocated to gain and shape quantization and 
define N g = 2 Ba and N s = 2 Ba . Further, let C g and C s represent the gain and shape codebook 
respectively: 



C 9 



C 9l J C 92 1 ' ' ' ) c 9n s 



(10) 

[c si , C S2 , • • • , C SJV __ J . (11) 

The product codebook can therefore be represented as, 

C = C g x C s . (12) 
The quantized gain and shape variables are computed as: 

g = arg min(j - c g ) 2 (13) 

CgGCg 

s = arg min I Is — cJ I 2 . (14) 

c s ec s 

The Lloyd-Max algorithm is the optimal solution to find the codebook for the gain of the 
channel vector with the MSE objective [|20l . We use the K-means approach, as described in [|2TT|. 
for numerical implementation of the Lloyd-Max algorithm. The optimal codebook of unit norm 
vectors with a Euclidean measure is not yet known. Therefore, we adopt random VQ to find the 
shape codebook. With this approach, the unit norm quantized shape vectors are randomly and 
independently distributed on the complex unit hyper-sphere in C M . 

The remaining question is, given B, what is the optimal choice of B s and B g l 

III. Distortion analysis and optimal bit allocation solution 
A. Design Objective 

Our main problem is to optimize the shape-gain bit allocation as formulated below, 

arg min £ [| Iz-S§l I 2 ] (15) 

Bs,Bg 

subject to : B s + B g = B , B s > , B g > , g e C g , s G C s . 
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Hamkins et al. [14] have shown that, for high resolution quantization (large B s and B g ), the 
distortion measure takes the following form [fT4ll : 

E[\\z-gs\\ 2 ]^E[(g- g) 2 ] + E [g 2 ] E [| |s - s| | 2 ] (16) 
^D g + E[g 2 } D 8 , (17) 

where, E [g 2 ] denotes the variance of the gain and D g = E [(g — g) 2 } is the gain quantization 
distortion. On the other hand, D s = E [| |s — s| | 2 ] represents the distortion due to unit-norm shape 
quantization. Since D g and D s are independent of each other in (fT7l) . the optimal bit allocation 
problem can be solved using the following three steps: 

1) Find D g , gain distortion, for a given B g . 

2) Find D s , shape distortion, for a given B s . 

3) Provide optimal bit allocation to minimize the overall distortion, i.e., E [g 2 ] D s + D g . 



B. Distortion due to Gain Quantization 

The distortion due to quantizing the gain is given by 

D g = E [(g - g) 2 ] = (r - g(r)f f g (r)dr. (18) 
Jo 

Here, g(r) is the quantized value of r and f g (r) is the probability density function (pdf) of the 
gain. Using Bennett's integral ( ifTTll . page-186), the distortion in ([TBI takes the form, 



where, N„ = 2 Bg and 



D 9 = ^\\f 9 (r)\\ k , (19) 
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ll/*(»0lli= / \fM*dr) . (20) 



Lemma 1: For Rayleigh fading and based on the pdf of the dominant eigenvalues of Wishart 
matrix and Jacobian transformation in [|22l and [|23l , we have, 

n /f(r )ii. = a^ F 3 (m±i\ , (21) 



4(L(e) - 1)!" 

where, L(e) = (M — e)(N k — e), M represents the total number of transmit antennas at the BS, 
N k denotes the number of receiver antennas of the k th user, e denotes the index of the ordered 
eigenvalues where represents the most dominant one, 1 denotes the 2nd most dominant one 
and so on. Finally, (3 = X e /L(e) where A e is the mean of the e th eigenvalue. 
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Proof: See Appendix A. ■ 

Using (fl9l ) and (I2TT) . the gain distortion at high resolution can be expressed as, 

A, = I ^jjll/,(r)||. (22) 
1 ' 3** 



16iV| (L(e) - 1)!' 
= K g 2- 2B % (24) 

where,if s = ^ (l^-i)^ 3 ( L ^3 +1 ) * s a constant with respect to B g . Equation (|24l) suggests that 
the gain distortion due to quantization is proportional to 2~ 2Bg . 

Figure [3] shows the distortion due to gain quantization of the dominant singular value of a 
2x2 MIMO channel. As the figure verifies, the analytical expression converges to the simulation 
result as B g increases. 

C. Shape Quantization Distortion 

This section focuses on the shape quantization distortion of a unit-norm vector in C M , in 
terms of the Euclidean distance. The Euclidean distance of two points in a C M plane has a 
one-to-one relation with the distance of two points in a IR 2M plane. Therefore, we can focus on 
quantization of unit- norm vectors in IR 2M instead of C M . 

Figure |4] shows a two dimensional view of the problem where OB = s, OA = s. Here, 
IMh = ||s||2 = 1. The Euclidean distance between s and § is defined by, d — ||s — s|| 2 . Define 
U.2M as the unit hypersphere in M. 2M . The surface area of U 2 m is given by |[T5l 

S(U 2M ) = 2MC 2M , (25) 

where, 

° 2M ~ r(M + i)' (26) 

Define the spherical cap V, i.e., the region ABC around s in Fig. SI as: 

V= (sGZY 2 m|||s-s|| 2 < d), (27) 

and let ZAOB = 9 be the angular distance between s and s. Since | \OA\ \ 2 = \ |s| \ 2 = 1, we have 
AD = sin(6 l ) and OD = cos(#). Also, since ||0-B|| 2 = ||s|| 2 = 1, we have BD = 1 - cos(6 l ). 
Therefore, 

AB 2 = AD 2 + BD 2 = sm 2 (6) + (1 - cos(^)) 2 = 2-2 cos(0). (28) 
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Here, if we define 6 = d 2 , we will have: 

6 = cos' 1 (1 - 0.56) . (29) 
The surface area of V is given by lfT5l . 

r 9 

S (V) = (2M - l)C 2M -i / sm 2M " 2 d<f>. (30) 



Now, if we assume a small spherical cap of radius d centered on s, the quantized vector can lie 
anywhere on this cap. Hence, 



Pr[||s-s|| 2 <6] = -|^-. (31) 
Using d25]), ©, and ® in (ED we get 

ivDi-»ii' < 6i = < 2M - jr 1 ™**"-** 

Ul " " J 2MC 2M 
Since all the quantized vectors are randomly chosen, the probabilities that the square of the 
Euclidean distance between any vector in the codebook and the corresponding channel is higher 
than 6, are independent of each other. Therefore, 



pr • 1 1 L (2M-1)C 2M -1 R 

Pr[ mm ||s — s»|| > b\ = 1 - 



cos" 1 (1-0.56) . 2M-2 x X 

sin 



(33) 



<6ti,JV.] \ 2MC 2M 

Hence, expected value of the distortion error due to shape quantization can be calculated as 
follows: 

£7(6)=/ Pr[min||s - s,-|| 2 > b]db. (34) 
Jo ie7V 

The limits of integration in (l34l) follow from the fact that the square of the Euclidean distance 
between two points on a unit radius sphere has a range of to 4. 
Lemma 2: 

£7(6) < £T s 2^r, (35) 

where, 

rr-l 

27T M T (2M_1 + ^ 

is a constant. 

Proof: See Appendix B. ■ 



^ = „ » fT , mTi , 7n • (36) 
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Figure \5\ shows that the upper bound of the shape distortion in (1351 ) has a fixed gap with 
respect to the simulation result. Therefore, we can safely approximate the shape distortion with 
the analytical expression in (|37l ). Thus, 



/,, ,,9\ / 7T 2 T(M) \ -2B S 

D s = E( s-s 2 ) « — - ; . 2**=* = K 8 2*n=* . (37) 

vii ii > 1 27r r (^f=± + i)y 

Z). Optimal Bit Allocation 

Having analyzed the individual terms in (TT71) , we are now able to answer the core question of 
this paper: the allocation of bits across gain and shape. In (fTTT) . the overall distortion measure 
was shown to take the following form, 

D = E [\ |z - gs\ | 2 ] w D g + E [g 2 ] D s . (38) 

Using the gain and shape distortion measures of (1241) and (1371 , D can be approximated as, 

D ss E [c/ 2 ] iT s 2-^T + ir 9 2" 2B « (39) 
w i? s 2"^fe + K g 2- 2{B - Bs \ (40) 

where i? s = -ft's-E [(? 2 ]. With these relations in hand, the optimal shape-gain bit allocation can 
be formulated as follows, 

B* = argmini? s 2-5ife + K g T 2{B - Bs) (41) 

Bs 

Theorem 1: 

The optimal bit allocation problem has the following solution: 

2M- 1 2M - 1 , / i?, 

S H log 9 ; : 

2M AM B2 \K g (2M-l 



Bs = ^v^B + _-_log 2 ( 7r77 ^ - ) (42) 



1 2M - 1 / K. 



B = B - B s = B - — log 9 I — r I . (43) 

9 s 2M 4M S2 \K g {2M - 1) ) V ' 

Here, K s and K g are the terms defined in the previous subsections. 
Proof: See Appendix C. ■ 

Note that, K s and K g in (l42j) and (|43T) depend on M but not 5. Therefore, as B goes to 
infinity, 

B s w 2M ~ 1 B (44) 
2M v ; 



1 

2M" 



5, ~ (45) 
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The analytical expressions of (1441) and ( 1431) can be intuitively explained as follows: The norm 
of a C M vector varies across a one dimensional line. However, the shape of a C M vector is 
uniformly distributed in the surface of a (2M — 1) dimensional hypersphere. Therefore, given 
2M number of bits to quantize a C M vector, one should expend approximately 1 and (2M — 1) 
bits to quantize the gain and shape of the vector respectively. It is worth noting that, from 
a different point of view and using a very different analysis, the work in lfT5l . [fT6l leads to a 
similar expression and explanation. However, this similarity is only for a high available feedback 
rate. 

Finally, the overall quantization error for a fixed feedback overhead takes the following form: 

D =K S 2~^ + K g 2~ 2Ba (46) 

( 2M-1 p , 2M-1 , / K a \\ of 1 p 2M-1 , / K s \\ 

-K 2 2M - 1 v 2M 4M g2 v if 9( 2M - 1 )// 4- JC 2 v 2M 4M ^V^gC 2 ^- 1 )// 



1 2M-1 



l0g * l^3T) J l^ 2 ^ - K g 2-^r) (47) 
=^2-^, (48) 
where, D c = log 2 ^ (2A/-i) ) {k s 2^^m — K" g 2~ 2 w i j is a constant. 

E. Overall Linear Precoding Algorithm 

In the previous section we derived the optimal allocation of available bits across gain and 
shape. Here we use this information to develop the overall linear precoding algorithm. The 
material exploits previous work in [|T9ll , 11241 . The algorithm steps are: 

1) A gain codebook of B 9 bits is generated based on the dominant singular values of a random 
Gaussian matrix and using the fC-means algorithm ETI . A shape codebook of 2 Bs random 
unit-norm vectors, uniformly and independently distributed in C M , is also generated. Both 
these codebooks are generated off-line and the codebooks are shared between the BS and 
the users. 

2) The BS sends common pilot symbols so that the receivers can estimate H fc . 

3) The receivers calculate the dominant singular values A k and the corresponding singular 
vectors V fc and form F fe = H fe V fe A fe . 

4) The receivers use the codebooks to quantize the gain and shape of the column vectors in 
Ffc and feedback the quantization indices to the BS. 
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5) 



The BS calculates the optimum virtual uplink power allocation matrix as, 

Q opt = min (a 2 + a ^ a A tr(J- 1 ), such that tr(Q) < P m 




max- 



Here, a\ = D as in (l48l) and J is calculated according to ©. 
6) The precoding matrix of the k th user is calculated as, TJ™ mse = J~ 1 Ffc V /qfc. Here, = 
[Qki, " " jQk L ] contains the virtual uplink power variables of the L k streams of the k th 



7) Using a recent result ||25l , the downlink transmit power variables are determined as, p = q. 
Here, q = ■ • • , qi] is the virtual uplink transmit powers of the L streams. 



This section provides the results of simulations to study the effect of shape-gain quantization on 
the performance of the MU MIMO linear precoding scheme described. We assume the following 
scenario in the simulation setup: the base station has two transmit antennas and serves two 
receivers in the downlink. Each receiver has 2 receive antennas and receives 1 data stream. The 
feedback overhead per user, B, is 16 bits. We show performance curves for different shape 
quantization bits, B s , in Figs. [6l [7] and [8] The corresponding number of gain quantization bits 
is given by, B g = B — B s . 

Figure [6] illustrates the effect of bit allocation on the quantization error and suggests that 
B s = 13 and B g = 3 are the optimal bit allocations for this scenario. The analytical results in 
(|42|) and (1431) lead to, B g = 2.6, B s = 13.4 which matches the numerical result. 

Fig. H plots the SMSE of the same system with the transmitter using 16-QAM. The figure 
shows that B s = 12 leads to the minimum SMSE. The SMSE performance of B s = 13, i.e., 
the optimal solution obtained from analytical results, is very close to that of B s = 12 bits. The 
minor difference between the simulation and analytical result stems from the fact that the derived 
gain distortion holds only for large number of gain quantization bits. Note that, B s = 16, i.e., 
quantizing the shape exclusively, leads to much higher SMSE. Therefore, optimal bit allocation 
across gain and shape feedback provides better performance in terms of SMSE. 

Fig. [8] shows that the bit allocation B s = 12 or B s = 13 also lead to better performance in 
terms of BER. If one uses all the bits for direction quantization, i.e., B s = 16, the effect of 
multiuser interference on the norm of the received signal cannot be rectified. This leads to the 
inferior performance of B s = 16. 



user. 



IV. Numerical Simulations 
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We did not plot the performance of all possible bit allocations, e.g., B s = to B s = 8, 
so that the figures look clearer. However, the performance trend of B s = 11 to B s = 9 bit 
clearly suggests that full gain quantization (B s — 0, B g — 16) will also perform much inferior 
to the optimal bit allocation in shape-gain quantization. Thus, optimal shape-gain quantization 
can improve over full gain or full shape quantization and lead to a lower BER in multiuser 
MIMO systems. In wireless ethernet ll26l systems, where a small number of bit errors may lead 
to the whole packet drop ll27l . optimal bit allocation in shape-gain quantization can significantly 
reduce the packet loss rate and save packet re-transmission time. 

V. Conclusion 

This paper studies the optimal bit allocation across gain and shape quantization in a MU 
MIMO downlink system by minimizing the SMSE of the system for a fixed feedback overhead 
per user. We show that the distortion due to gain and shape quantization are proportional to 2~ 2Ba 

2B 3 

and 2 2M1 respectively, suggesting that, in the asymptotic region of high feedback overhead, the 
number of shape quantization bits should be approximately (2M — 1) times than the number of 
gain quantization bits. The analysis and importance of bit allocation is borne out by the simulation 
results that show significant worse performance for the usual approach (in MU MIMO downlink 
systems) of only quantizing the gain or shape but not both. 

Our work with respect to the gain distortion calculation is quite general, since the gain 
quantization distortion of other distributions like Rician, Nakagami and Weibull fading can 
also be calculated using Bennett's integral. However, the optimal bit allocation results might 
be different from the Rayleigh fading case considered in this paper. 

Appendix A 
Proof of Lemma 1 

The authors of |[22l| have provided the following pdf of the eigenvalues of a MIMO channel, 



Here, A e denotes the e eigenvalue of the Wishart matrix (i.e., H H or HH ). e denotes the 
index of the ordered eigenvalues. L(e) = (M — e)(N k — e). (3 is a constant whose value is given 




(49) 
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through the following equation, 

— T7~y (50) 
L{e) 

Here A e is the mean of the eigenvalue. (l49l) provides the probability distribution function of the 
eigenvalue of the Wishart matrix, A e . In our proposed algorithm, we are trying to quantize g, 
the singular values of the Gaussian matrix H. Now, A e = g 2 . 

Using Jacobian transformation E3l . the probability distribution of the singular values of the 
Gaussian matrix can be found as follows, 

i ( r 2\L(e)-l / 2 

Therefore. 



-» r > = (L(e) - !), ffLW ^{-j) 2r - < 5 » 



HAM"! = (I(iFT^^> (T^"* (-Q *)'■ <52) 



Using standard mathematical tables of Il28l ( P - 380, eqn - 662), we find 



r oo r (a±i 

x n exp (—ax p ) dx = — , +1 ( . (53) 
3 pa^~- 



Comparing (|53l) with (|52l , we find, n = 2L ^ 1 , a = p = 2. Therefore, 



3 3£ 

2L(e)-l 



2L(e)-l 



r 3 exp ( -— jdrj= 2L ■ - (54) 



2 (4 



°° 2L(e)-l 



r= ^ rap (-5) *•) = \ m ^ r (^ir 1 ) <55> 

'~ r ^ exp = I (3 «M«» +lr 3 (41+1) . (56) 

Using (HU) in d52"b, we get, 



3 x 3 L ^(3 3 / L(e) + 1 
4(L(e) -1)! V 3 



(58) 
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Appendix B 
Proof of Lemma 2 



Using (|33]>, 



Pr [mm | |s - Sjjf > 6] = 1 - — 



cos 1 (1-0. 5b) . 2M-2 

sin . 



A? 



-.v \ 2MCW 

iV 



(59) 



/cos" 1 (1-0.56) 
sin 2iW -odo ) (60) 

/cos-^l-O.Sfe) \ 
<f) 2M - 2 d(j)J (61) 

= (l - K 2 (cos" 1 (1 - 0.56)) 2A/_1 ) * . (62) 

In (l60l) . we assumed = ( - 2A ^^"~ 1 . (|6TI) follows from the fact that, given a large number 
of quantization vectors, i.e., at high bit rate, the complementary cumulative distribution function 
(CCDF) is significant only for smaller values of 0. For these smaller angles, we can assume 
sin</> « (p. Equation (l62l) follows from assuming K 2 = 2 m-i ' 

Figure |9] compares the simulated shape distortion with the original and approximate analytical 
shape distortion of a 2 x 1 C M vector. Here, the CCDF of the original and approximate analytical 
expressions are superimposed with the simulated CCDF. Hence, (|60~1) and (I6TT) accurately model 
the actual distortion. This justifies the transition from (l60l) to (|6"TT) . 

Now, using (l34l) . we find, 

E{b) = [ Pr[mm\\s - Si\\ 2 > b]db (63) 
Jo ie7V 

(l - K 2 (cos' 1 (1 - 0.5b))™- 1 ) d& (64) 

= 2 / (1 - K 2 B 2M ' l ) N sin{e)dB (65) 
</o 

^2 [\l - K 2 6 2M - l ) N 0d6 (66) 
io 

^2 f 1 (1 - K 2 9 2M - l ) N 9d9 (67) 
io 

= 2 jf ( E (f) (-l)^0 i(2Af " 1)+1 J ^ (68) 
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= 2 ^Mhy~2- (69) 

The transition from (1631 to (|64l) can be explained as follows: the similarity between (1601 and 
(|6T| ) holds only for smaller values of b since sin <p ^ 4> for larger <p. Therefore, although the 
square of the Euclidean distance between two random unit norm vectors can vary from to 
4, (|62l holds only for a smaller range of b. At the presence of a large number of codewords, 
the squared distance between the original and the quantized channel takes large values with a 
negligibly small probability. Therefore, we can truncate the range of b as long as the CCDF of 
the original function is negligible outside the range, i.e., the limited range of b does not have any 
significant affect on the calculation of the expected value of the distortion. Using this analysis, 
in (|64|) . we use a as the truncated range, i.e., we assume that b can vary from to a. 

In (|65l) we assumed, 9 = (cos -1 (1 — 0.56)). Therefore, ip = (cos -1 (1 — 0.5a)). Since only 
smaller angles of 6 contribute to E(b), we assumed sin6> m 6 in (l66l) . In (|68l) . we assumed ip = 1 
to simplify the other calculations. 

Fig. [10] justifies the approximations that we used in the derivations of shape distortion cal- 
culation. Here, approxl and approx2 denote sin(6>) ~ 9 (ref: eq. |66]) and rj 1 (ref: eq. 1671) 
respectively. As Fig. [10] shows, the three curves are superimposed with each other. Therefore, 
our justifications are valid for high bit rate quantization. 

Applying (fj = ( ~ l)i( ~ N)i , where (-JV)< = ^tSt ^3' ^ takes the following form, 



r(-A0 

X M-1)'(--V),(-D'A^ 2 ^ HjOgjg 

4^ z!(i(2M-l) + 2) 2 M- 1^1(2 + ^) 

2 JV! 



A -l)*(-jV) t (-l)*ir 

, n i!(i(2M-l) + 2) - 2 _ u _, 



K^- 1 (71) 



^3 (72) 



r(iv + i + dbi) 
ivr(iv)r (jf±i) 
r(iv + |S±i) 



K 3 (73) 



IB was found using ( ED, 6.6.8). In (|72l> . we assumed if 3 = K 2 2M1 . d74]) was obtained using 

»r(fc) 

T(a+fe) 



the relation between the gamma and beta function, f3(a,b) = ^+1? iBTTl. Following a similar 
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work in 0, we find, 



/ 2M + i\ R r(2 B )r(i + ?x f T ) 

AT/3 N, — = 2 B ) J (75) 



= r ( 2 * + 1 > (77) 
The preceding inequality in (l76l) is justified by the following reasoning: due to the convexity 



of the gamma function [3] and the fact that = T(2) = 1, T(x) < 1 for 1 < x < 2 . Let, 

V = 2S + 2aI=i' * = 1 - 2jJ=I' so that ' 2/ + 1 = 2 s + 1, y + 1 = 2 B + 1 + ^fn- By applying 
Kershaw's inequality for the gamma function ll32ll . 



r(y + t) / , t 



Using d78l), 



ny+1) <(V + i) 1 V«>0,0< ( <1. (78) 

r(2 B + i) / R 2 i \^ 

-2 

^ \ 2M-1 

-2B 

<2^i. (81) 



Using (1811) and the value of K 3 we find, 

N 



^ z!(i(2M-l) + 2) \2MC 2M ) ' K ) 

Using the values of C2M-1 and Cim one can obtain, 

E(b) < K s 2^, (83) 



( 2M-1 \ 2A/-1 

2 ^ r |2M-i +1 j J is a constant with respect to B s . 

Appendix C 
Proof of Theorem 1 

Taking the 1st and 2nd order derivatives of (l40l) . we find, 
dD 



dB s 
d 2 D 

d? Be 



R.Qn2)2-&* (-^ry) + K g Qn2) (2" 2 ^)) 2 
^(ln 2) 2 2"^r (-^- T ) 2 + K 9 (2 In 2) 2 ( 2 - 2 ^) 



(84) 
(85) 
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From (|85l ), J^p > 0. Therefore, the optimal bit allocation problem is convex Il33i Now, equating 
the 1st derivative to be zero, 

k " -2^ = K g 2- 2 ( B - B ^ (86) 



2M 



2 -2B+2B s + ^ til (87) 

K g (2M-l) 1 } 

2MB, „ 1, / K s \ 

B +n^g 2 (——^-—) (88) 



2M - 1 2 toi \K g (2M - 1 



2M - 1 2M - 1 . / K 



= B + log 9 — ; — r . (89) 

2M AM * 2 \K g (2M - 1) ' 



Therefore, at the optimal point, 



2M - 1 2M - 1 , / K 

E> _i_ 1 . . (y I ; 

2M 4M &2 \K g {2M-l) 



B s = -^v^B + _ rT _log 2 ( - ) (90) 



1 2M-1 ( K s \ 
B *=2M B - - W l0g2 U(2M-1) J ' (91) 
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Fig. 2. Gain-shape product quantization 




Fig. 4. Shape quantization block diagram 
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-a — Simulated Distortion 

1 — Analytical Expression of the upper bound 
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Fig. 5. Comparison of the simulated distortion with the theoretical upper bound (2x1 complex vector) 
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Fig. 6. Effect of bit allocation in the quantization of the product of dominant eigenvalue & the corresponding eigenvector of 
a 2 x 2 MIMO channel 
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Fig. 7. Effect of bit allocation in the SMSE of 16-QAM system, M = 2, N = [2 2], L = [1 1], B = 16 



M = 4, N = [2 2], L = [1 1], 16 bit per user, 16-QAM 
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Fig. 8. Effect of bit allocation in the BER of 16-QAM systems, M = 2, N = [2 2], L = [1 1], B = 16 
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Fig. 9. Comparison of the original and approximated complementary cumulative distribution function of the shape distortion 
of a 2x1 vector (10 bit quantization) 
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Fig. 10. Justification of the approximations used in Shape distortion calculation 



